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ABSTRACT 

In this study reacting proton with intermediate and heavy nucleus as a target ( 64 Ni, 67 Zn, 77 As, llh Rh , 1 11 Cel, 114 Cd, 
I67 Ho, I 69 Tm, l86 W) the experimental data of cross sections have been published in Exfor library as a function of proton 
energy. We have calculated the cross section of the above mentioned data and results have been obtained by using 
(Matlab-8.3 2014a) program. The stopping powers have been calculated from Zeigler formula by using SRIM-2013 with 
the results of cross sections to calculate the neutron-yield for reactions, and also comparing between cross section for each 
one of those reactions with published experimental data of cross sections at Exfor library as a function of proton energy to 
show compatibility, as well as those reactions to be used in the production of radioisotopes such as ( 64 Cu , 67 Ga, 75 Se, '" 'Pd, 
UI In, II4m In, l65 Er, l69 Yb, 176 Re ), We also did found comparisons between neutron yields of the mentioned reactions to 
choose the best reactions of highest neutron yield. 

KEYWORDS: Cross Section (Excitation Function), Stopping Power, Neutron Yield, Data Evaluation, Radioactive 
Isotopes 

INTRODUCTION 

The first major advance of particle accelerators occurred in 1934 with the invention of the cyclotron, it is machine 
used to accelerate charged particles to very high speeds, and it became possible to create the nuclear instability that is a 
prerequisite for radioactivity. By directing a beam of fast moving charged particles at a nucleus target, it induced a reaction 
that results in the formation of a radioisotope with a short half-life [1]. Particle accelerators and, in particular, cyclotrons, 
were very important in the preparation of radioisotopes during the years from 1935 to the end of World War II. After that 
reactors were used to produce radioactive elements, but the use of accelerators for this purpose became less common. As 
the techniques for using radiotracers became more sophisticated, it was clear later that reactor produced radioisotopes 
could not completely satisfy the growing demand and, therefore accelerators were needed to produce new radioisotopes 
that could be used in new ways in both industry and medicine [1]. The study of nuclear reactions of proton with 
intermediate and heavy nucleus as a target which has a high neutron yield due to the advantages of the ( p,n ) reactions to 
produce radioisotopes as ( 64 Cu , 67 Ga, 75 Se, I03 Pd, 111 In, II4m In, 165 Er, 169 Yb, 186 Re). They are very important because it 
provides a base for a wide range of technical applications, particularly using the reactions induced by intermediate energy 
protons; it is possible to produce directly radioisotopes which are used in medicine and industry, recent decades, where 
widespread uses of diagnostic and therapeutic radioisotopes had occurred. Depending on the type of radiation, the 
diagnostic isotopes are classified into two groups: /? + , ff emitters ( 64 Cu, "'"Re, etc.) used in Positron Emission Tomography 
(PET), E.C and y-emitters ( 67 Ga , 77 Se, '"'Pd, 111 In, 114m In, I65 Er, >69 Y,) used in Single Photon Emission Computed 
Tomography (SPECT). The total cross section of such a production is also important in accelerator technology from the 
point of view of radiation protection safety [2]. 
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NUCLEAR REACTION 

The g-value 0 f the X(a,b)Y reaction is defined as the difference between the final and initial kinetic energies and 
is given by [3,4.5]: 

Q = [M a +M x -(M b + M y )]x 931.5 (1) 



The threshold energy for a nuclear reaction is defined as the smallest value of a projectile’s kinetic energy at 
which a nuclear reaction can take place [5, 6], 
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NUCLEAR CROSS SECTION 



To study nuclear reactions, it is necessary to have a quantitative measure of the probability of a given nuclear 
reaction. In which the cross section of nuclides is the effective area presented [7]. When an accelerated charged particle 
interacts with a target nucleus a nuclear reaction takes place, ultimately leading to a stable or radioactive product nucleus. 
A nuclear reaction is characterized by a cross-section, describing the probability that a particle can interacts [6, 8], 
The cross section has the units of area as it’s the square of nuclear radius. The standard unit for measuring a nuclear cross 
section (a) is the barn (b), which is equal to (lbarn=lCT 24 cm 2 ) or (lbarn=10 -28 m 2 ). The reaction cross section data 
provides information of fundamental importance in the study of nuclear systems. The cross section is defined by [9]: 

rr = *- (3) 

/ 

Where R is the number of reactions per unit time per nucleus. I is the number of incident particles per unit time 
per unit area. In general, a given bombarding particle and target can react in a variety of ways producing a variety of light 
reaction products per unit time. The total cross section is than defined as [5]: 

Vunal = (4) 

Where a , is the partial cross section for the process. 

STOPPING POWER 



The energy of charged particle loses per unit path length of the material it traverses. Generally speaking, any 
charged particle can have either electronic, nuclear, or gravitational interaction with the particles of the material through 
which it passes. However the gravitational interaction is too low to be of any significance and is generally ignored. The 
total stopping power is then just the sum of the stopping powers due to electronic and nuclear interactions [10, 11, 12]. 
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Where S e i ec tronic is the electronic stopping power, S nuc i ear is the nuclear stopping power. The electronic stopping 
power is always much larger than the nuclear stopping power. For most practical purposes the nuclear component of the 
stopping power can also be ignored as it is generally only a fraction of the total stopping power. For particles such as 
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electrons, this statement is always valid since they are not affected at all by the strong nuclear force. 

Electronic Stopping 

Electronic Stopping Power (S e ) the beam particles hitting a target get slowed down by interactions 
with the electrons until it they are in thermal equilibrium with their surroundings. As a consequence one gets a wide 
neutron spectrum if the projectiles are completely stopped in the target, even if all nuclear reactions were two-body 
reactions. The loss of kinetic energy in a nuclear encounter would be much larger [4], In the scope of this work, the 
electronic stopping powers were calculated using the Ziegler formulae [11]. Reliable data are available for many elements 
over a wide range of energies. However, in order to obtain values for all elements over a continuous range of proton 
energies, the stopping power as a function of proton energy E (keV) and the target atomic number Z 2 . S e is assumed to be 
proportional to E 0 ' 45 for E < 25 keV, except for Z 2 < 6 where it is proportional to E 0 25 , for 25 keV < E < 10 MeV 
[12,13]. 



(— )-(- — ) + (- ) 
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S Low (Low energy stopping) is 

5 Low = A 1 £^+A 3 £^ (7) 



Sffigh (High energy stopping) is 

A, ln(^- + A v £) 
c - E 



Where E(KeV) energy of proton and the coefficients A, for each Z 2 , available from SRIM [12,13]. For 10 MeV < 
E < 2 GeV. 
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For use in the parameterization Proton stopping, S t ( 4 H), in various target materials for 10 kev < E < 100 MeV, 
High-energy proton stopping Se ~St , Proton stopping power for 1 kev < E < 10 GeV and Z 2 < 92, based on a combination 
of theoretical calculations and experimental data (for 20 keV < E < 1 MeV) [14,15]. 



Nuclear Stopping Power (S n ) 

A beam of charged particles bombarding the neutral atoms of a target interacts with the atomic nuclei and atomic 
electrons of the target. The ratio of the energy lost in interaction with the atomic electrons, to the energy lost in the 
interaction with the atomic nuclei (2m p / m e = 4xl0 3 ) [10]. 

Thus the energy lost by interaction with the nuclei is negligible compared with that lost by interaction with the 
electrons. For solid targets with a given thickness there will be the more nuclear interactions and hence the electronic 
stopping is comparable with the nuclear cross section [11]. The theory of ion-solid interactions is well established and it 
has been discussed in detail. Where the authors also described the computer program SRIM. Which can be used for 
calculating stopping power and range for ion target combinations with Z<92 [12, 14, 16]. For sufficiently high projectile 
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energies, S,~ S e . S n , in units of eV/10 15 atoms/cm 2 , for any projectile with energy E (KeV), is given by [12,14]. 



SJE) = 



8.462Z 1 Z 2 M 1 5„(g) 

(M x +M 2 )(Zj 0 ' 23 + Z 2 ' 23 ) 



( 10 ) 



Where the reduced ion energy, s , is defined as 

32.53 M l M 2 {E!M l ) 



Reduced Ion Energy = £ ■ 



Z l Z 2 (M l +M 2 )(Z,°- 23 +Z 2 °' 23 ) 



( 11 ) 



Where M, and M 2 are the projectile and target masses (amu), and Z, and Z 2 are the projectile and target atomic 
numbers. For e < 30 keV. 
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For s > 30 keV , unscreened nuclear stopping is used, and S„ (s) simplifies to 
Inf 



Sj£) = 



2e 



(13) 



NEUTRON YIELDS 



The neutron yield (Y n ) detected per incident particle (proton, deuteron alpha), for an ideal thin, and uniform target 
and mono-energetic particles beam of incident energy E b is given by [17]. 



K = (N d x)<J (E b )n(E b ) 



(14) 



Where N d is the real number density of target atoms, x is the target thickness, o is the reaction cross section and r| 
is neutron detection efficiency. If the target is sufficiently thick, for a target which is not infinitesimally thin, the beam 
loses energy as it passes through the target, if the target is sufficiently thick, and there exist one atom per each molecule 
(i.e.,/=l) and taking (E')=l, then the resulting yield is called the thick-target yield which is given by [17,18,19]. 

v - {e ‘'= n }M? e <15) 

Where E lhr is the reaction threshold energy, o(E) is the cross section, dE/dX is the incident particle initial energy, 
N is the atomic number of target per unit volume, which is defined as follows [20]: 

N = W P N « (16) 

A 



Where, w is the abundant in the combination, is the combination density, A is the mass number, N A is the 
Avogadro’s number. Thus by measuring the yield at two closely spaced energies (El) and (E2), one can determine the 
average value of the integrand over this energy interval as follows [20]: 



y{e 2 )—y{E\)= 
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Where ( E b ) is the average of (E,) and (E2). If a (E) are available in the literature as a function of projectile energy 
(E b ) for natural elements, then the neutron yield can be calculated using eq. (13). If neutron yield is available as a function 
of projectile energy (E b ), then eq. (13) can be used to calculate a (E) as a function of (E b ). Thus consequently one can 
calculate the neutron yield by using eq. (17), for natural elements and if only one stable isotope is available in nature [21]. 

RESULTS AND DISCUSSION 

These reactions are the important interactions that lead to the radioactive isotopes production such as (Copper-64, 
Gallium-67, Selenium-75, Palladium-103, Indium-111, Indium-1 14 m , Erbium-165, Ytterbium-169 and Rhenium-186). 

The jg AT (p , n) “ Cu Reaction 

The cross sections of the ^ Ni (p ,n) Cli reaction have been published as a function of proton energy by 

Adam [22], Avila [23], Szelecsenyi [24], Levkovskij [25] in Exfor library. We’ve used the energy range (1.0-24.5MeV) 
lower than the threshold energy (2.4956MeV) within the step of (0.5 MeV), the cross section is directly proportional with 
proton energy until it reaches to (11.0 MeV), it’s the maximum value of the cross section (o Ma x=742.95 mb), then it 
decreases with the increase of proton energy as shown in figure 1 . 

The f Q Zn(p ,n) 3j Ga Reaction 

The cross sections of the Zn(p ,n) b lGa reaction have been published as a function of proton energy by 

Szelecsenyi [26], Hermanne [27], Levkovskij [28] and little [29] in Exfor library. We’ve used the energy range (4- 
29.5MeV) higher than the threshold energy (1.8105MeV) within the step of (0.5 MeV), the cross section is directly 
proportional with proton energy until it reaches to (11.0 MeV), it’s the maximum value of the cross section (a Max =739.09 
mb), then it decreases with the increase of proton energy as shown in figure 2. 

The 33 As(p,n) USe Reaction 

The cross sections of the 33 As(p,ri) reaction have been published as a function of proton energy by 

Levkovskij [30], Brodovitch [31], Johnson [32], and Albert [33] in Exfor library. We’ve used the energy range (2.0- 
59.95MeV) higher than the threshold energy (1.6693MeV) within the step of (0.8MeV), the cross section is directly 
proportional with proton energy until it reaches to (9.2MeV), its maximum value of the cross section (GMax=692 mb), then 
it decreases with the increase of proton energy as shown in figure 3. 

The "yRh (p , ri) '^Pd Reaction 

The cross sections of the “ 3 Rh (p , n) Pd reaction have been published as a function of proton energy by Sudar 

[34], Hermanne [35], Harper [36] and Johnson [37] in Exfor library. We’ve used the energy range (2.305-39.055 MeV) 
higher than the threshold energy (1.3385MeV) within the step of (0.75 MeV), the cross section is directly proportional with 
proton energy until it reaches to (9.805MeV), it’s the maximum value of the cross section (o Ma x=509.93 mb), then it 
decreases with the increase of proton energy as shown in figure 4. 

The "g 1 Cd(p,n) ' 4 ' 9 '/« Reaction 
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The cross sections of the Cd ( p , n ) "Jjn reaction have been published as a function of proton energy by 

Tarkanyi [38], Marten [39], Skakun [40] in Exfor library. We’ve used the energy range (3.8-43.8MeV) higher than the 
threshold energy (1.6596MeV) within the step of (0.8MeV), the cross section is directly proportional with proton energy 
until it reaches to (11.8MeV), it’s the maximum value of the cross section (cMax=755mb ), then it decreases with the 
increase of proton energy as shown in figure 5. 

The "g Cd(p,n) U4 ^In Reaction 

The cross sections of the Cd ( p , n) lH ^In reaction have been published as a function of proton energy by Said 

[41], Tarkanyi [42], Zaitseva [43] and Wing [44] in Exfor library. We’ve used the energy range (4.4-62.4MeV) higher than 
the threshold energy (2.2486MeV) within the step of (l.OMeV), the cross section is directly proportional with proton 
energy until it reaches to (11.4MeV), it’s the maximum value of the cross section (oMax=284.81 mb), then it decreases with 
the increase of proton energy as shown in figure 6. 

The '^Ho(p,n) ‘“£r Reaction 

The cross sections of the l ^Ho(p ,n) ^ Er reaction have been published as a function of proton energy by 

Tarkanyi [45], Beyer [46] in Exfor library. We’ve used the energy range (5.5-37.5MeV) higher than the threshold energy 
(1.1674MeV) within the step of (0.6MeV), the cross section is directly proportional with proton energy until it reaches to 
(10.9MeV), it’s the maximum value of the cross section (oMax=150.1mb), then it decreases with the increase of proton 
energy as shown in figure 7. 

The u f 9 Tm(p,ri) 'f () Yb Reaction 

The cross sections of the l6 ? 9 Tm(p ,ri) 1 %Yb reaction have been published as a function of proton energy by 

Tarkanyi [47], Sonnabend [48], Spahn [49] in Exfor library. We’ve used the energie range (3.3-44.7MeV) higher than the 
threshold energy (1.6910MeV) within the step of (0.9MeV), the cross section is directly proportional with proton energy 
until it reaches to (11.4MeV), it’s the maximum value of the cross section (G\ tlix = 291 mb), then it decreases with the 
increase of proton energy as shown in figure 8. 

The 18 7 ®W(/?,?t) Rc Reaction 

The cross sections of the ( p , n ) l8 x Re reaction have been published as a function of proton energy by 

Hussain [50], Xiaodong [51] in Exfor library. We’ve used the energy range (5.0-17.8MeV) higher than the threshold 
energy (1.3700MeV) within the step of (0.4MeV), the cross section is directly proportional with proton energy until it 
reaches to (9.0MeV) and it’s the maximum value of the cross section (o Max = 39.231 mb), then it decreases with the 
increase of proton energy as shown in figure 9. 
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FigUrel Figure 2 




Incident proton Energy (MeV) 



Figure 3 




Figure 4 





Figure 5 



Figure 6 



www.iaset.us 



editor @iaset. us 



84 



Shafik Y. Qastin & Huda M. Tawfeeq 






Figure 9 



CONCLUSIONS 

From the results of cross sections and neutron yield for (p,n) reactions with many of intermediate and heavy 
nucleus, we have observed in this study many conclusions for this study as per below points: 

From the results of cross sections for (p,n) reaction with many intermediate and heavy nucleus targets to produce 
radioisotopes with short half life, those radioisotopes are producing by cyclotron more demanded [52], 

Copper-64: has a half-life (T 1/2 =12.7 h) and decay properties (]3 + , 0.653 MeV [17.8%] (3‘, 0.579 MeV [38.4%]) 
[54]. These characteristics make it useful for both high resolution PET imaging and targeted end radiotherapy. In addition, 
its electron capture decay associated with Auger emission gives more efficient cell killing when this radioisotope is 
deposited in the cell [53]. The well-established applications in the molecular imaging of cancer [54]. 

Gallium-67: has a half-life of (T 1/2 =3.217 d) decays to stable 67 Zn by electron capture. Its decay emissions 
include gamma rays of 93.3 keV (37.0%), 184.6 keV (20.4%) and 300.2 keV (16.6%). Gallium behaves in the body in a 
similar way to ferric iron. It is commonly used as a trivalent citrate compound for nuclear medicine imaging, and is a 
valuable agent in the detection and localization of certain neoplasms and inflammatory lesions. [52] This radioisotope 67 Ga 
is well known and widely used in the field of nuclear medicine. 67 Ga has become one of the most frequently employed 
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cyclotron produced radioisotope over the last two decades and is a widely used single photon marker for detecting the 
presence of malignancy and for diagnosis of inflammatory diseases. [55] y-emitters from 67 Ga used in Single Photon 
Emission Computed Tomography (SPECT) [55]. 

Selenium-75: has a half-life of (Ti / 2 =1 19.78 d) decays to stable 75 As by electron capture 100 %. Its decay 
emissions include gamma rays of 120.8(17.2%), 136(58.5%), 265(58.6%), 279.5(24.9%) [57]. used in the form of seleno- 
methionine to study the production of digestive enzymes [59]. 77 Se radioisotopes are widely used in medicine, industry and 
agriculture. 75 Se is used in high activity brachytherapy, assessment of pancreatic exocrine function, study of bile acids and 
evaluation of illeal function, industrial radiography and as a tracer in the assessment of chemical, biochemical, biophysical 
processes, metabolic researches and agricultural studies [57]. 

Palladium-103: has a half-life of (Ti / 2 =16.9991 d) decays to stable 103 Rh by electron capture 100 %. Which de- 
excites by means of a heavily converted internal transition. As a result of both processes (EC and IT), X-rays and Auger 
electrons are emitted which are ideally suited for brachytherapy (its electron capture decay resulting in abundant emission 
of Auger electrons and low energy X-rays (20-22 keV)) [98, 58]. '" I’d used to make sealed seeds (brachytherapy) 
permanent implant seeds for early stage prostate cancer [59]. 103 Pd is extensively used in the treatment of prostate cancer 
and ocular melanoma, and is mostly applied in the form of sealed seeds (brachytherapy) [58]. 

Indium-111: has a half-life of (T, /2 =2.83d) decays to stable UI Cd by electron capture 100 %. There are two 
prominent gamma rays, one at 171.3 keV and one at 245.4 keV [52]. Indium-111 octreotide, a radiolabeled somatostatin 
analogue, binds to somatostatin receptors, which are very common in several cancers. Indium-111 may, therefore, be 
useful for the visualization of metastases in cancer patients. [52]. 11 'in Used for specialist diagnostic studies, brain studies, 
infection and colon transit studies [59]. 

Indium-114 m : is nucleus of the nucleus of 1 1 'In in the first energy level excited. ' l4 '"ln has a half-life of 
(T i / 2 =49.5 1 d) decays to stable by electron capture (E.C) with positron ( + fi) 3.25 % and decay to ll4 In by Internal Transition 
(IT) 96.75%, II4m In used in radio-immunotherapy, therapeutic radioisotope. The use of positron emitters radionuclides such 
as ( 64 Cu, ll4m In, 123 I, I49 Tb, 195m Hg, etc), for radionuclide therapy is questionable because the interaction of positrons in 
tissues results in the 511 keV annihilation photons. These constitute a major contribution to the absorbed dose in the 
surrounding healthy tissues. Nevertheless, positron emitters may be valuable for the selection of radio-immunotherapy 
candidates by confirming tumor targeting and/or estimating radiation doses [60]. 

Erbium-165: has a suitable half-life of (Ti / 2 =10.36 h) decays to stable '""Ho by electron capture 100 %. Erbium- 
165 is one from very few radionuclides exists that decay exclusively by EC-mode without any accompanying radiation, 
16S Er is one of them. Auger electrons are emitted by isotopes that decay by electron capture (EC) or have internal 
conversion (IC) in their decay of these isotopes a cascade of very low energy electrons is emitted [61]. 

Ytterbium-169: has a half-life of (T 1/2 =32.0 1 8 d) decays to stable 169 Tm by electron capture 100 %. I69 Yb can be 
used for cerebrospinal fluid studies in the brain [59]. Many methods to produced of 169 Yb in nuclear reactor or by 
accelerators cyclotron. The average cross section of the 16S Yb(n,gamma) 169 Yb process is very high as compared to that of 
the 169 Tm(p,n) I69 Yb reaction. However, the abundance of l(lS Yb in natural ytterbium is only 0.13 %. In practice. In cyclotron 
production via the I69 Tm(p,n) I69 Yb reaction, the cross section is lowest then reaction in reactor , but the abundance of I69 Tm 
in natural 100% [62]. 
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Rhenium-186: has a half-life of (T 1/2 =3.7183 d) decays to stable /,S6 Wby electron capture (E.C) 7.47 % and decay 
to 186 Os by electron emitted ('/?) 92.53 % with E Milx p =1.1 MeV, Ey =137 keV [63]. 186 Re Used for pain relief in bone 
cancer, beta emitter with weak gamma for imaging [59]. 186 Re is a beta and gamma emitter which is used in metabolic 
radiotherapy and radio immunotherapy (RIT). After labelling with hydroxethylidene diphosphonate (HEDP) is a new 
compound used for the palliation of painful skeletal metastases. Make it a useful radioisotope for the treatment of small 
tumors in the body. One of the most common applications of I86 Re is palliation of painful bone metastasis caused by 
prostate or breast cancer, and it is also useful for treating painful arthritis. [64] In addition, the y-ray emitted with E y =137 
keV (branching ratio of 9.5%), which can be imaged by Single Photon Emission Computed Tomography (SPECT) to 
determine the distribution of the radioactivity in the body [63]. 

In figure 10, we’ve noted abundant reactions with high neutron yield value when the proton is interacting with 
intermediate nuclei target such as ( 64 Ni, 67 Zn , 75 As, W3 Rh, 1,1 Cd, U4 Cd ) compared with heavy nuclei target ( I65 Ho , I69 Tm, 
m W). 

In figure 11, the " As(p,n ) USe reaction has high neutron yield compared with 16 lTm(p ,n) >69 0 Yb reaction, as the 

stopping power of incident proton on nucleus target 77 As is lower than I69 Tm and also the 73 As(p,ri) ^Se reaction has a 

high cross section as shown in figure (3) with dependence on equation (17), when incident proton on intermediate mass 
nuclei target 75 As has a high neutron yield compared with heavy mass nuclei I69 Tm target. 

When incident charged particle such a proton on intermediate and heavy mass nuclei 103 Rh and "’ Ho targets, We’ve 
noted the Rh[p,n ) ^Pd reaction is high neutron yield compared with l87 Ho(p,n) reaction, as shown in figure 12 
because w ’Rh target has low stopping power of Incident proton compared with "' Ho target, as well as the abundance of 
llh Rh and " Ho is 100%, the ]< "Rh{p,n) '^Pd reaction has a high cross section as shown in figure (4) with dependence on 
equation (17). 




Figure 10: Neutron Yield of Intermediate and Heavy 
Nuclei as a Target 




Figure 11: Neutron Yield of Intermediate and Heavy 
Nuclei as a Target 
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Figure 12: Neutron Yield of Intermediate and Heavy Nuclei as a 

Target 
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